function[fday,hourly_ts, ts_meanremoved] = wdc_read_ts(fid,element),
[data,data1,fday1] = wdc_read(fid,element);

hourly_ts = reshape(data',1,size(data,1)*size(data,2));

st=1;
en=24;
for i = 1:size(data,1),
    fday(st:en) = fday1(i)+(0.5:23.5)./24;
    st = st+24;
    en = en+24;
end;
